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Observations of the temperature structure function 

2 

parameter, C T , over the ocean 



K. L. Davidson, T. M. Houlihan, C. W. Fairall and 

G. E. Schacher 



Environmental Physics Group 
Naval Postgraduate School 
Monterey, CA 93940 



Abstract 

Observational results from shipboard measurements 

of temperature fluctuations, and mean wind, temperature, 

and humidity, are compared with existing expression for 

2 

the surface flux and height dependence of C T . Surface 
flux estimates are obtained from bulk aerodynamic for- 
mula. Temperature fluctuation data are selected to 

minimize a salt-contamination effect which causes in- 

2 

creases in ' emperature variance. Predictions for C T 
based on surface flux scaling agree within 20%, except 
for near neutral and large unstable conditions where the 
disagreement can be attributed to measurement problems. 



1 . 



Introduction 



The atmospheric marine surface layer is that region in the 
first 50 meters or so over the sea which is characterized by 
large turbulent vertical transfer of momentum, heat, and mois- 
ture. Progress has been made in specifications of surface wind 
fields and hydrostatic stability for this region since these 
are parameters used in estimating boundary fluxes in numerical 
atmospheric prediction and for predicting changes in the upper 
part of the ocean. Another interest in this region is the 
quantification of turbulent effects on optical wave propagation. 
The optical effects occur at scales which are smaller than 
those associated with boundary fluxes. This involves specifica- 
tions of fluctuations in velocity, temperature, and humidity 
for scales approaching those responsible for viscous dissipa- 
tion of turbulent kinetic energy and molecular diffusion of 
scalar (temperature and humidity) quantities. 

Existing measurement capabilities and turbulence theory 
enable observations of the small scale properties of interest, 
but the complexities and cost of such direct observations in- 
hibit them on routine bases. Therefore, there is a need to 
examine relationships between the small scale properties and 
the more feasibly measured parameters: the mean winds, tem- 
perature, and humidity. 

Descriptions of the small scale turbulent properties have 
not been as complete nor in the quantity for the overwater 
regime as for the overland regime (described by Champagne, et 
al . , 1977),. The overwater regime requires separate considera- 
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tion because of evidence of influence on the airflow by the 
surface waves (Davidson, 1974) and the importance of the water 
vapor content of the air. The water vapor effects are mani- 
fested both in atmospheric stability and in optical propaga- 
tion . 

For optical waves, the largest contributors to the refrac- 
tive-index fluctuations are temperature fluctuations (T'). 
Friehe, et al . (1975) have shown that another significant 
contributor may be the temperature-humidity fluctuation co- 

variance. They observed a contribution of 24% to the refrac- 

2 

tive index structure function parameter, C n , from humidity 
effects during a time when the temperature-humidity covariance 
was positive. 

The intensity of temperature fluctuations for scales from 

a few millimeters to several meters can be characterized by 

2 

the temperature structure function parameter, , defined as 



T(x+r) are temperature at two points on a line oriented normal 
to the mean wind direction separated by the distance r. 

The temperature structure function parameter can also be 
related to small scale processes of turbulent kinetic energy 
dissipation and heat diffusion by the following expression 




( 1 ) 



D T (r) is defined as D T ( r ) = [T(x) - T(x+r)] 2 where T(x) and 



C, 



2 



n e 



-1/3 



X 



( 2 ) 



T 



2 



where n is a constant equal to 3.2, e is the dissipation rate 
of turbulent kinetic energy and x is the rate at which the 
temperature inhomogeneities are diffused or "smeared” out. 
Equation 2 arises from two parallel expressions for the one- 
dimensional variance spectrum for temperature fluctuations in 
the inertial subrange. 

S T (k) = a C T 2 k~ 5/3 (3) 

and 

S T (k) = 6 e -1/3 x k' 5/3 (4) 



where a and 3 are constants equal to .25 and .8 respectively. 

The latter expression was formulated by Corrsin (1951). 

This paper describes results from a comparison of a very 

2 

substantial amount of C T data obtained over the open ocean with 
theWyngaard et al. (1971) predictions based on bulk surface 
layer measurements. Friehe (1976) performed a similar compari- 
son in which fewer data were used and in which a different 

bulk stability parameter was used. He concluded that reason- 

2 

able prediction of C T could be made from bulk formulations. 

2 

2 . Surface Laver Bulk Aerodynamic Relations for C^ , 

2 

Wyngaard, et al. (1971) related C T to surface fluxes with 
Eq. 2 using separate expressions for e and x based on overland 
measurements. The general expression was 



P 2 „2/3 

T - 2 ~ = 5 - l 



(5a) 
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where 



f(€) = 



5 < 0 



(5b) 

(5c) 



4.9(1 - 71 ) 



-2/3 



4.9(1 +2.4 5 2 ^ 3 ) 



S > 0 



It is noted that these interpolation formulae still lack sub- 
stantial verification for the overwater regime. 

These formulae agree with the expected asymptotic scaling 

n 

(Wyngaard, 1973) for C T under very unstable (local free 
convection) and very stable (Z-less) stratifications. The 
appropriate scaling formulations in the extreme stability limits 
are 



free convection: 



where 



Z-less : 



C T 2 z 2 / 3 



= constant 



( 6 ) 



T f = (TQ 2 /gZ ) 1/3 ; Q q = 



w' 0 



C T 2 L 2 / 3 



= constant 



(7) 



2 

Therefore, the asymptotic predictions are for C T to vary as 
Z _4y/3 for free convection and to be independent of height 
(Z-less) for very stable conditions. 

In this study, boundary flux parameters in Eq. 5 , (T* and 
L) were calculated using bulk aerodynamic formulae. Bulk 
aerodynamic formulae are those which relate boundary scaling 
parameters to the wind speed at a level and the temperature 
and humidity difference between that level and the surface. 
Their use for indirect estimates of the boundary fluxes has 
received considerable attention in overwater observational 
experiments because of difficulties in making eddy flux, and 
profile measurements. 
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The bulk aerodynamic formulae at a height Z for the veloc- 



ity and potential temperature scaling parameters are: 



u * = V u 


(8) 


1 


t* = c 0 2 (G z - e s ) 


(9) 



where S refers to the surface value and Cq and C 0 are drag co- 
efficients. The drag coefficients depend upon height, 
stability, and wind velocity. Similarly, for the sensible 
heat flux 

-vPlF = C R U (0 Z - 0g) ' (10) 

where 

C H - < C D C 0 )4 (11) 

Since several publications give detailed descriptions of 

the bulk method with various experimental values for the drag 

coefficients, we shall restrict this discussion to a few 

points that must be considered when relating overwater measured 
2 

values of C^, to values of T + and £ based on bulk method 
calculations. In particular, we must note that T* is directly 
dependent upon the value of C 0 2 (Eq. 9). Using conventions 
developed by Liu (1978), the neutral stability temperature drag 
coefficient is given by (in this convention , (a,pK Z/T*) (30/3Z) 
is the dimensionless profile function) 

C Q / = a T K[£n(Z/Z oT )]~ 1 (12) 
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where k is von Karman ' s constant, the ratio of heat and 

momentum transport dif fusivities (at neutral stability) and 
Z ,p is the appropriate roughness length for the temperature 
profile. Since Wyngaard, et al.'s expression for f(£) was 
based upon direct measurement of T* from the heat flux and 
Reynold's stress, their curve was not dependent upon a choice 
of Cg^. In view of this, we will select the value of C SN that 
yields the best fit to our data, treating this as an indirect 

measurement of CL... 

0N 

The dependence of overwater neutral momentum drag coeffi- 
cients , , on wind speed has been the objective of many in- 

vestigations. For our analyses, we chose a recently suggested 
representation by Kondo (1975), Table 1. 

The stability dependence of and is defined on the 
basis of the following general expressions 



c = ^DN 

D " [1 - CO - 1 C DN * f ^)] 2 

c = C 0N^ 

[1 -(a T K) _1 C 0N ^ 4> 2 (£)] 2 



(13) 



(14) 



where <1^(5) and are stability corrections to wind 

and scalar profiles. Formulations for ^(O and ip 2 (£) have 
been presented by Businger (1973) and ours were similar except 
for the use of a T normalization (Liu, 1978). 

The dependence of the drag coefficients on stability leads 
to computational considerations in estimating £ since the co- 
efficients which define the latter are themselves functions 



6 



of £ as indicated by the following bulk expression 




(1 - < 1 C “DN ^(O) 2 
(l - ( a rpK ) ^ C“Q^ ^2^^ 



( 15 ) 



where 



gZ C9 N 2 



o 



j- A0 + . 18AQj 



(16) 



All parameters are in the MKS units with Q in gm/kg. 

5 is the equivalent of a first guess of £ based on bulk 

differences and neutral values for the drag coefficients. For 
each value of wind velocity, we calculated C^, allowing cal- 
culation of The final value of 5 was found by solving 

Eq. 15, iteratively. An example set of solutions is shown in 
Fig. 1 for U = 7.6 m/sec. 

Friehe (1976) effectively defined £ as being equal to 

2 

when he compared overwater values with bulk parameter 

predictions based on Eq. 5. He used a constant C DN value and 
two different values of C^, depending on the value of UA0 . 

The differences between the solid curve and dashed lines in Fig. 
1 represent the difference which would arise in neglecting stab- 
ility effects in the specification of £. 

3. Experimental Arrangements 

The data were obtained during several experiments conducted 
during a four year period (1974-1978), over an extensive ocean 
area and aboard three different ships (listed in Table 2). The 
instrumentation was designed for profile measurements of both 
mean and turbulent parameters. In all experiments, these were 
made at multi-levels on two masts separated spatially on the 
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forward part of the ships (Fig. 2). A comprehensive descrip- 
tion of the system is the subject of another paper (Houlihan, 
et al. , 1978). 

Mean temperature was measured with quartz oscil- 
lator thermometers, accurate to .01°C. Mean humidity was mea- 
sured with Li Cl sensors, accurate to about 3%RH. 

Wind velocity was measured with cup anemometers, accurate to 

o 

about 5%. Temperature fluctuations (for C T ) were measured 
with ac Wheatstone bridges (3.0 KHz carrier), using paired 
resistance probes with 2 . 5ym diameter platinum wires. The 
resistance probes were mounted on wind vanes to maintain correct 
alignment with the relative wind direction. 

Two features of the measurement arrangement significant in 
the interpretation of results are the aspiration of mean tempera- 
ture sensors and the surface temperature measurements. The 
quartz temperature sensors and humidity sensors were contained 
within aspirated weather shelters. The shelters had two purposes: 
to protect the sensors from the marine environment and to reduce 
radiation effects. Radiation effects were of special concern, 
since some of the sensors were above the water and some were 
above the deck of the ships, which are heated by the sun and 
radiate strongly. Because of radiation from the deck, we modi- 
fied the shelters, placing a shield at the base of the unit to 
reduce radiation from below. Even with these precautions, system- 
atic temperature corrections for the aspirators (during day and 
night operations) have been determined from calibrations with a 
glass dewar type aspirator. The daytime correction .6C and the 
nighttime correction is .4C (subtracted from measured temperature). 
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The sea surface temperature sensor (a quartz probe) 
was attached to a plastic hose so that it floated within 
one inch of the surface. On all ships, its location was 
aft of the midship. The location on all of these ships, 
underway and holding position, was definitely influenced 
by the ship's wake. 

The mean data were continuously averaged and periodic- 
ally logged by a microprocessor developed at the Naval 
Postgraduate School and adopted for the R/V ACANIA by Plunkett 
(1977). Mean values were defined for 15 to 30 minute periods. 
Most periods were 30 minutes long, the planned duration, but 
shorter periods occurred due to such events as change in the 
ship's course. 

4. Data Evaluation and Error Anal ysis 

2 

Comparison of C T results and available expressions are 
made with T* and L values estimated from 10 meter wind, 
temperature, and humidity and surface temperature values. 
Initial plans were also to estimate T^ and L values from 
mean profile data, as evident from previously described 
instru nentat ion arrangements. The latter results are not 
presented because of inconsistencies which have been attribu- 
ted to difficulties in defining small scalar and wind gradi- 
ents. Although these same errors existed in the mean data 
used for bulk estimates, the accuracy requirements of the 
bulk method are much less stringent. 
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2 

The 10 meter mean, as well as C T values, were obtained by 

interpolating from levels above and below, if available, or by 

extrapolating on the basis of expected height variations, if 

only one level was available. With these specifications, 

and £ were computed for the 10 meter level, using Eqs. 9 and 15. 

Two aspects of the measurements require further discussion 

for this study. These are the influence of salt contamination 

on resistance wire sensors, and uncertainties in the surface 

temperature specification. Although the latter may not be sig^ 

nificant for some applications, it is significant in our results, 

because the air-surface temperature difference occurs as a 

2 2 

squared quantity in scaling C T by T* . 

The influence of salt contamination of resistance wires 

and thermocouples, and hence, all temperature fluctuation 

measurements in the marine environment, was recently described 

by Schmitt, et al. (1978). These effects were discussed by 

2 

Friehe (1977) in his comparison between oyerwater values 
and bulk aerodynamic predictions and led him to limit the data 
considered. Schmitt, et al. (1978) suggested that erroneous 
temperature fluctuations occur due to water vapor absorption 
and evaporation on the salt nuclei attached to the wire. This 
occurs in conjunction with fluctuations in ambient HH. The 
effect is manifested by the occurrence of spikes in the tempera- 
ture trace. These spikes can cause overwater temperature 
variance spectra to exhibit non -Kolmogorov slopes. 
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In an effort to eliminate the salt effect, the resistance 

wires were washed or changed at frequent intervals during the 

latest experiment (CEWCOM-78). Mean ratios of immediately 

before and after washing versus relative humidity appear in 

o 

Fig. 3. We observed that C T values from contaminated wires 
are larger for an RH range from 50% to 85%, only. Presumably, 
the effect is diminished at RH above 85% because humidity fluc- 
tuations are less intense due to the small RH gradient. Their 
effect, as described, is within the system accuracy. Presumably, 
the effect is diminished at RH below 50% because the salt nuclei 
attached to the wires are not activated, hence, humidity fluc- 
tuations do not cause significant absorption or evaporation. 

We tested these interpretations by computing the ratio of 
2 2 

measured C^, to computed C T (using Eq. 5 and bulk parameter 
estimates of T* and £) with data from all experiments except 
CEWCOM-78. Mean values of these ratios versus RH appear in 
Fig. 4. The individual points are mean values for 5% RH 
intervals and the error bars are expected errors of the aver- 
ages (± a//N , N = number of points) within each interval. 

The number at the top of the error bar is N. These data exhibit 
values near 1 for RH values above 85%. There are slight in- 
creases at 95% and at 85% and there is a significant increase 
at 80%. The ratios are larger at low RH values, compared to 
corresponding result from clean wires in CEWCOM-78. 

Differences between the results in Fig. 3 (washing and re- 
placement) and the results in Fig. 4 (from the other experiments) 
are the larger values of the ratio at very low RH values and the 
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lower values above 95%. The former result is not unexpected, 
because wires exposed for longer periods without cleaning as 
they were for the other experiments would be sufficiently con- 
taminated to respond to large humidity fluctuations, even if 
they were not activated. The latter result could also be 
attributed to the degree of contamination and the tendency of 
salt particles to grow into water droplets at high RH. Con- 
taminated wires at large RH values experience a loss of fre- 
quency response, resulting in a lower sensitivity to tempera- 
ture fluctuations. 

Based on the above results, we decided to consider only 

2 2 
those C,p data measured at RH values above 82.5%. results 

obtained with lower RH values were excluded except from periods 

immediately following sensor replacement. All data from 

CEWCOM-78 which were from clean sensors were considered. 

2 

We used intermediate C T and A0 results to examine un- 
certainties in the sea surface temperature. These results 

2 

appear in Fig. 5 where C T data were selected on the basis of 
RH considerations and A0 values have been adjusted for system- 
atic aspirator corrections. If the measured surface tempera- 
tures were the "skin" temperatures, we expect the minimum in 
2 

C,p to occur near A0 equal to zero. This is the case for both 

night and day. However, we estimate that our A0 values are 

subject to uncertainties as large as .5°C, based on the dis- 

2 

placement and broadness of the daytime C^, minimum. This is 
primarily due to an inability to ensure that the sea surface 
sensor actually measured the skin temperature. Liu (1978) 
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suggests how the latter problem could have been treated with 
interface flux equilibrium considerations. However, his 
suggested procedures were not applied to these data because 
they apply to a bucket temperature which was not measured by 
our sensor. 

The results of an analysis of measurement errors appear in 
Table 3. These are the errors introduced by instrument accura- 
cies, noise and uncertainties in certain calibration values. 

2 

In the case of C^, , most of the uncertainty is due to the salt 
contamination effect. In the case of Tg, the sensor itself 
is accurate to .01°C, but one cannot be certain that the 
sensor is in fact measuring the surface temperature. 

5 . Results 

Dimensionless temperature structure parameter (DTSP) versus 
£ results appear in Fig. 6, in the format used by Wyngaard, et 
al. (1971). Again, error bars are the uncertainties of the 
means and the numbers are the number of points defining the 
mean. For this comparison, C Q ^ was determined by assuming that 
Z _ was independent of velocity and had a value corresponding 

_3 

to Z q for a typical drag coefficient, Cp N = 1.3 x 10 . The 

2 

appropriate value of Cg^ for a measured C T was estimated from 
the following equation. 



'0N 



= c 



- < 

c 2 z 2 / 3 _ (£n Z/Z oT ) -1 t(; 2 (C)]‘ 



TUT 



-) 



A0‘ 



(17) 



The value of C Q . T selected was that from the CEWCOM-78 data, 
u JN 



since these were least influenced by salt contamination. 



Although it was necessary to assume an initial value of 
Cq^ to obtain the values of £ , the final results (Fig. 7) 
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are fairly insensitive to the initial choice. A minimum |£| 

value of .05 was used in this determination in view of the 

uncertainties listed in Table 3. The value of Cg^ selected 

and, hence, used to compute T* and £ for the results in Fig. 6 
-3 A 

was 1.3 x 10 . Using C RN = (C^ Cq^) 2 , we compare this C 0N 

value with some obtained from measurements of (Table 4), 

-3 

assuming a typical C D ^ value of 1 . 3 x 10 if not given. 

DTSP results are also presented with logarithmic £ axes, 
Fig. 8. This format enables examination of the distributions 
of the large number of results near £ equal to zero in a more 
illustrative perspective. Also, results at £ values less 
than -7 which were not included in Fig. 6, can now be viewed. 

This format demonstrates the exceptional agreement between 
the results and the prediction over the range, -.03 < £ < -7. 
Also, disagreement between observed and predicted results can 
be readily related to previously mentioned measurement uncer- 
tainties. Disagreements at the minimum |£| intervals, for 
both stable and unstable sides, are attributed to errors in 
specifying A0 as well as inherent noise in the C ^ measurement 

and recording systems. Disagreement at large |£| values on 

q 

the unstable side is attributed to increased C T values due 
to the salt contamination effect, since these data happen to 
have corresponding relative humidity values near 85%. 

DTSP results based on asymptotic scaling for unstable 
(Eq. 6, free convection) and stable (Eq. 7, Z-less) conditions 
appear in Figs. 9a and 9b, respectively. Predictions based on 
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Eq. 5 (a and b) appear as solid curves. For the unstable side, 
the prediction was approximated as 

c T 2 Z 



which approaches a value of 2.7 for large |?|. Eq. 18 is not 
correct when q*, as well as T* , determines £. When q* is 
negative (typical marine conditions), the observed values 
should be less than predicted by this expression. This is 
the case for these results. For the stable side, the predic- 
tion is 

C 2 L 2 / 2 

— — = 4.9 [e’ 2/3 + 2.4] (19) 

T* 2 

which approaches a value of 11.8 for large ?. 

In summary, the DTSP results agree very well with the 
predictions under unstable conditions and reasonably well 
under stable conditions. Under stable conditions, they are 
slightly less than the predictions. The determined Cq N value 
agrees reasonably well with those reported. DTSP 
disagreement under neutral conditions (|?| < .05) is attribu- 
ted to uncertainties in measuring the air-sea temperature 
difference and inherent system and measurement noise. The 
latter would include the salt contamination effect, since 
any humidity fluctuation would cause apparent temperature 
fluctuations. Disagreement under very unstable conditions, 

£ < -7, was due to the fact that these data occurred with 



2/3 



= 9.9 



2/3 



l-7£ 



(18) 
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RH values, = 85% , not totally immune to the salt contamina- 
tion effect. 

The actual significance of the slight disagreement under 
stable conditions is somewhat uncertain. However, it could 
also be attributed to the salt contamination effect. During 
stable conditions over the ocean, humidity and temperature 
fluctuations are negatively correlated. Hence, dry-warm or 
moist-cold eddies are causing the fluctuations. In the case 
of dry-warm parcels, contaminated sensors would be experienc- 
ing a cooling due to evaporation, but a compensating heating 
due to the temperature change would also occur. In the case 
of cold-moist parcels, the reverse would occur. 

DTSP results based on asymptotic scaling appeared to 
approach the predicted constants at the limits. Agreement 
for these results was expected, on the basis of agreement in 
other DTSP results, because Monin-Obukhov scaling was used 
in determining T^. This arose when Q q , in Eq. 6, was deter- 
mined from T + and U* values which were estimated on the basis 

of Monin-Obukhov scaling. 

2 

The multi-level C T data enabled an examination of 
predicted height variation, under different stability condi- 
tions. These results appear in Fig. 10 where A is the 

coefficient obtained for the best fit linear relation between 
2 

observed Log C T values and Log Z values, 

Log C T 2 = A Log Z + B. (20) 



16 



The curves are predictions obtained by determining dLog 

2 

C(p /dLog Z from Eq. 5 and approach limiting values of -4/3 
and 0 on the unstable and stable sides and a value of -2/3 
at the neutral point. 

Fewer periods were available for this examination than 

2 

for the preceding C T results. This was because probes at 

upper levels, which were not readily accessible, were not 

always replaced when they became inoperative. Also, we 
2 

measured C,p at a single level during CEWCOM-78 because 

of frequent probe washing and replacement. 

Results in Fig. 10 indicate that height dependence of 
2 

observed values agreed quite well with the prediction. 

This is particularly the case for the stability variation 

of the dependence. A majority of the mean values indicate 
2 

that C T decreased somewhat less rapidly with height than 
the prediction. This result could be attributed to system 
noise which contributed the same at all levels and would reduce 
the apparent relative difference (decrease) between levels. 
Evidence for this is the larger (relative) -disagreement at near 
neutral (unstable) intervals, where the signal to noise ratio 
is the smallest. 

6 . Conclusions 

2 

Comparisons of a substantial amount of overwater C^ data 

have been made with predictions of the existing empirical ex- 

2 

pressions. Measured overwater C T data, which were selected 
to minimize the salt contamination influence are described 
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very well by Wyngaard, et al's (1971) expressions, even 
with bulk surface flux estimates. An estimated valu 

_3 

of 1.3 x 10 agrees with lower values obtained by others. 

Disagreements exceeding measurement accuracies (20-30%) occur 

2 

for near neutral conditions where C,p is insignificant in most 
applications, and for large unstable conditions. These dis- 
agreements are attributed to measurement difficulties rather 
than to the validity of the existing expressions. Agreement 
between observed and predicted height distributions over the 
full stability range provide further verification of the exist- 
ing expression. 

In conclusion, a typical (but not the best) comparison of 

2 

observed and predicted C,p values during a 24 hour period in 
CEWCOM-78 is presented in Fig. 11a to further demonstrate the 
overwater applicability of the expressions. The predicted 
values based on bulk estimates of surface fluxes are within 
the accuracy of the measurements in describing the observed 
values. Coincident bulk |?[ estimates for the period appear 
in Fig. lib. 
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Table 1 



C D n versus wind speed (10 meter) 
from Hondo (1975). 



U ms 



-1 



C DN X 10 



.3 


- 2.2 


1.08 


X 


u - 15 


2.2 


- 5.0 


. 77 


+ 


.086 x 


5.0 


- 8.0 


. 87 


+ 


.067 x 


o 

00 


-25.0 


1.2 


+ 


.025 x 



U 

U 

U 
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TABLE 2. Summary of Experiments 



EXPERIMENT 


SHIP 


LOCATION 


PERIOD 


NUMBER OF 
PERIODS 




R/V ACANIA 4 


Monterey Bay 


All Months 
1973-76 


98 


CEWCOM-76 1 


R/V ACANIA 


E. Pacific 


Sep-Oct 76 


285 


NRL - EOMET 2 


USNS HAYES 5 


N. Atlantic 
& Mediterranean 


May-Jun 77 


269 




USNS KANE 6 


Mid-Atlantic 


Feb-Mar 78 


38 


CEWCOM-78 3 


R/V ACANIA 


E. Pacific 


May 1978 


200 








Total 


890 



"''Cooperative Experiment West Coast Oceanography & Meteorology-1976 
o 

Naval Research Laboratories, Electro-Optical Meteorology Cruise 

3 

Cooperative Experiment We st Coast Oceanography & Meteorology"1978 

4 

Operated by Department of Oceanography, Naval Postgraduate School, 
Monterey, California 

5 

Operated by Naval Research Laboratories, Washington, D. C. 

6 

D Operated by Naval Oceanographic Office, Washington, D. C. 
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Table 3. 



Error analysis for measured , 
intermediate , and scaling parameters 



Measured 


T 

X S 


T 

Z 


U RH 


C ^ 


RMS Error (±) 


. 3°C 


. 2°C 


8% 3% 


30% 


Intermediate 




A9 


AQ 




Accuracy (1) 




. 5°C 


.4 gm/Kg 




Scaling 


u* 


T* 


Q* 


5 


Accuracy (±) 


10% 


o 

to 

o 

o 


.02 gm/Kg 


.03 or 30% 



24 



Table 4. 



Empirically determined sensible 
heat drag coefficients, 
x 10 3 



Sources 

Pond, et al. (1974) 
Dunckel, et al . (1974) 
Muller-Glewe and Hinzpeter 
Smith (1974) 

Friehe and Schmitt (1976) 
This study 

^( ) denotes assumed 

o 

Excludes high wind data 





C HN 


C DN 


n 

CD 




1.5 


1.45 


1.55 




1.5 


1.60 


1.40 


(1974) 


1.0 


( 1 . 30 ) 1 


.77 




1.2 


1.02 


1.41 




.91 2 


(1.30) 


.64 




1. 3 


1.30 


1. 34 


value ; 


other 


given by 


r authors 
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Figure 1. 


Ratio of E based on stability dependent drag 

coefficients ( E ) to E based on constant drag 

coefficients (E ) versus E . 

o o 


Figure 2. 


Schematics of shipboard mounting arrangements, 
(a) R/V ACANIA, (b) USNS HAYES, and 
(c) USNS KANE. 


Figure 3. 


2 

Ratio of C T (before to after replacement or 
washing) versus RH; (From CEWCOM-78). 


Figure 4 . 


2 2 

Ratio of measured C^, to computed C,p (from Eq. 5) 
versus RH, for all except CEWCOM-78 data. 


Figure 5. 


Crp"" versus A0 (a) daytime (0800-1900 LST), 
(b) nightime (1900-0800 LST). 


Figure 6. 


Mean DTSP versus E results. Solid curve is 
prediction based on Eq. 5. 


Figure 7. 


C Q .. versus E, results. C Qvr based on Eq. 17 and 
9h 0N * 

CEWCOM-78 results. 


Figure 8. 


Same as Figure 6 with logarithmic E axis. 


Figure 9. 


DTSP based on asymptotic scaling. (a) free convection 
curve is approximate prediction, Eq . 18 and, 

(b) Z-less, curve is prediction, Eq. 19. 


Figure 10. 


2 2 

Height variations of C T versus E • A(C^, ) is defined 

by Eq. 20, solid curves are prediction based on Eq. 5 


Figure 11. 


Results from 24 hour period during CEWCOM-78 

(a) Observed and Predicted C for 10 meters 

(b) Observed %. 1 
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10.0 
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Figure 1. Ratio of £ based on stability dependent drag 
coefficients (£) to 5 based on constant drag 
coefficients (£ ) versus £ . 




Figure 2. Schematics of shipboard mounting arrangements, 
(a) R/V ACANIA, (b) USNS HAYES, and 
(c) USNS KANE. 
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o 
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Figure 3. Ratio of C T (before to after replacement or 
washing) versus RII; (From CEWCOM-78). 



o 
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Figure 4. Ratio of measured C T to computed C T (from Eq. 5) 




2 

Figure 5. versus A0 (a) daytime (0800-1900 LST), 

(b) nightime (1900-0800 LST). 
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h ro 



100 




Figure 6. Mean DTSP versus £ results. Solid curve is 
prediction based on Eq. 5. 
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I 



Figure 7. Cg^ versus £ results. Cg^ based on Eq. 17 and 
CEWCOM-78 results. 
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Figure 8. Same as Figure 6 with logarithmic £ axis. 
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Figure 11. 



Results from 24 hour period during 

(a) Observed and Predicted C^, for 

(b) Observed £ . 



CEWCOM-78 
10 meters 
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